Height Fluctuations and Intermittency of V2O5 Films by Atomic Force Microscopy 
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The spatial scaling law and intermittency of the V2 O5 surface roughness by atomic force microscopy 
has been investigated. The intermittency of the height fluctuations has been checked by two different 
methods, first, by measuring scaling exponent of q-th moment of height-difference fluctuations i.e. 
C' q —< \h(xi) — h(x-2)\ q > and the second, by defining generating function Z(q,N) and generalized 
multi- fractal dimension D q . These methods predict that there is no intermittency in the height 
fluctuations. The observed roughness and dynamical exponents can be explained by the numerical 
simulation on the basis of forced Kuramoto-Sivashinsky equation. 
PACS: 52.75.Rx, 68.35.Ct. 
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Due to the technical importance and fundamental in- 
terest, a great deal of effort has been devoted to under- 
standing the mechanism of thin-film growth and the ki- 
netic roughening of growing surfaces in various growth 
techniques. Analytical and numerical treatments of sim- 
ple growth models suggest, quite generally, the height 
fluctuations have a self-similar character and their aver- 
age correlations exhibit a dynamic scaling form [1-6]. 

Vanadium pentoxide, V2O5, has been the subject of 
intense work because of its diverse applications in cat- 
alytic oxidation reactions, cathodic electrode in solid 
state micro-batteries, windows and electrochromic de- 
vices as well as gas sensors and to be of interest for 
transmittance modulation in smart windows with poten- 
tial application in the architecture and automotive. Also 
V2O5 is a low mobility semiconductor, and having pre- 
dominantly an n-type. Electrons are the charge carriers, 
and an increase in the carrier density in accompanied by 
reduction in oxygen concentration in the lattice [7-9] . 

This work aims to study the roughness and dynamical 
exponents ( \ an d z ) an d t ne intermittency of the V2O5 
films by the atomic force microscopy. Wc measure the 
height-difference moments C q {l = \x\ — 2:2 1 ) =< \h(x\) — 
h(x2)\ q > and show that they behave as ~ |a;i — X2|^ 9 . 
The obtained £ g is a linear function of q. We also in- 
troduce the generating function Z(q, N) and generalized 
multifractal dimension D q and show that D q behave also 
as a linear function of q. These observations indicate 
that the height fluctuations arc not intermittent, ft is 
also argued that the measured roughness and dynami- 
cal exponents belong to early-time scaling regime of the 
noisy Kuramoto-Sivashinsky equation. 



V2O5 layers were grown on the polished Si(100) sub- 
strate by resistive evaporation method in a high vacuum 
chamber. The pressure during evaporation was 10 -5 torr. 
The thickness of the growing films was measured in situ 
by a quartz crystal thickness monitor. We performed 
all deposition at room teperature, with deposition rate 
about 10 — Yhnmjmin. However, during the film deposi- 
tion the substrate temperature raised to T ~ 60C°, due 
to the radiation effect of alumina boat. The substrate 
temperature was determined using chromcl/alumcl ther- 
mocouple mounted in close proximity of samples. Surface 
composition of samples was measured by Auger electron 
spectroscopy (AES) using a 3 keV electron beam and 
a cylindrical mirror analyzer (Varian model 981-2607). 
The surface topography of the films was investigated us- 
ing Park Scientific Instruments model Autoprobe CP. 
The images are collected in a constant force mode and 
digitized into 256 x 256 pixels with scanning frequency 
of 0.6 Hz. The cantilever of 0.05 N m _1 spring constant 
with a commercial standard pyramidal S^A^ tips has 
been used. A variety of scans, each with size L were 
recorded at random locations on the V2O5 film surface. 
In order to determine the structure of the deposited V2O5 
films, we have performed XRD measurements of sam- 
ples. The spectra for the as-deposited V2O5 films grown 
showed that the V2 O5 thin films are amorphous. For 
thick films (d > 200nm) two very broad weak peaks 
were observed representing the growth of (001) and (002) 
peaks of the V2O5 orthorombic structure [7]. The yel- 
low color of deposited oxide films and their optical trans- 
mission spectra indicates the existence of V2O5 structure 
rather than other vanadium oxide phases [8]. AES anal- 
ysis of the V2O5 samples showed V and O peaks at the 
surface of the deposited films. The stoichiometry of the 
vanadium oxide films was calculated from the ratio of O 
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FIG. 1. AES spectra from V2O5 surface. The O/V ratio is 
2.5 ±0.1 (by considering the elemental sensitivity factors). 

to V Auger peak heights by considering the elemental 
sensitivity factors. The O/V ratio was 2.5 ±0.1 indicat- 
ing the formation of stoichiometric vanadium pentoxide 
at the surface of the thin film ( see Fig.l). With these 
observations we conclude the film composition is nearly 
stiochiomctric[9] . 

III. RESULTS 

Figure 2 shows the variation of the surface morphol- 
ogy for different growth times ( or thickness) of processed 
samples by AFM. As deposition proceeds, the size of 
mountains and valleys grows until the system approaches 
to a stationary state in which the thickness is about 
150mn. 

The quantitative information of the surface morphol- 
ogy can be derived by considering a sample of size L, 
and defining the mean height of growing film h and, its 
roughness w through [10]: 

1 f L/2 

h(L,t) = - dxh(x,t) (1) 

L J-L/2 

and 

w(L,t) = ({(h- hf)Y /2 - (2) 

(■ • •) denotes an averaging over different realization (sam- 
ples). Starting from a flat interface (one of the possi- 
ble initial conditions), it was conjectured by Family and 
Vicsck [11] that a scaling of space by factor b and of 
time by a factor b z (z is the dynamical scaling expo- 
nent), re-scales the roughness w by factor b x as follows: 
w(bL,b z t) = b x w(L,t), which implies that 

w(L,t)=L*f(t/L z ). (3) 

If for large t and fixed L (t/L z — > 00), w saturates then 
f(x) — > const., as x — > 00. However, for a fixed large 




FIG. 2. AFM surface images (all 0.5 x 0.5^m 2 ) of V2O5 
films with thickness of (a) 20nm , (b)100nm, (c) 260 nm. 
(from up to down) 



L and 1 << t « L z , one expects that correlations 
of the height fluctuations are set up only within a dis- 
tance t x / z and thus must be independent of L. This 
implies that for x « 1, f(x) ~ x^ with /3 = xl z - 
Thus dynamic scaling postulates that, w(L,t) ~ fr for 
1 << t « L z and ~ L x for t >> L z . The rough- 
ness exponent x and the dynamic exponent z charac- 
terize the self-afiine geometry of the surface and its dy- 
namics, respectively. We measure the exponent x from 
equal-time height-height correlation function defined as 
C 2 (l = \x\ - x 2 \) =< \h(xi) - h(x 2 )\ 2 >. Here /i(x) is 
the surface height at position x on the surface relative to 
the mean surface height. 

In order to determine the roughness exponent x — 
£2/2, we consider the roughness of the samples with 
thickness 20nm, lOOnm and 260nm. Fig. 3 shows that, 
for the thickness 20nm and lOOnm, the scaling behav- 
ior exists only for the scaling region ~ 4nm to 2Qnm, 
but for the stationary state sample, < h >= 260nm, 
there is a cross over for the scaling exponent x = £2/2 
in I = I* w 22nm. The exponent is £ 2 i f° r the length 
scales Anm to the scale and for the length scales 
25nm < I < bOnm, the fluctuation is determined by 
another exponent £22- The measured values for the ex- 
ponents £21 and 62 are 1.64±0.06 (xi = 0.83±0.03) and 
0.58 ± 0.08 ( X2 = 0.29 ± 0.04) respectively. The mea- 
sured roughness exponent x of the samples 20nm, lOOnm 
and 150nm are 0.71 ± 0.04, 0.77 ± 0.03 and 0.82 ± 0.03, 
respectively. As shown in Figs. 4 and 5, we note that 
the correlation length I* is about 22nm. Therefore, we 
are dealing with the statistical properties of the corre- 
lated scaling surfaces. As discussed in [22], such systems 
exhibit the sampling-induced hidden cycles (log-periodic 
fluctuations [23]). Such oscillatory behavior will dimin- 
ish when the sampling size is sufficiently large. The os- 
cilation amplitude approaches zero to within an order of 
5 = * /ML), where L and M are the sampling size 
and the number of independent curves which are aver- 
aged. We have determined each exponents by averaging 
on the eight AFM images. In our averaging it appears 
that there is no log-perodic property for the height-height 
correlation function. Therefore the surfaces of V2O5 are 
self-affinc. 

The existence of the cross over scale I* in the station- 
ary state is not observed only in the second moment C 2 . 
In Fig. 4 log-log plot of C3 vs I for the three samples are 
presented. It is evident that for the stationary state sam- 
ple, the scaling exponent of the third moment C3 has also 
a cross over in I*. The measurement shows that in the 
stationary state the exponents £3^ behave as £3,; = 3/2£ 2i 
with i = 1, 2. 

We also examine the scaling behaviour of the q-th mo- 
ment of the height-difference C q =< \h{x\) — h(x2)\ q > 
and show that all of the moments at least up to q = 20 
behave as \xi — x 2 \^ q ■ Fig. 5 shows £ g vs q for 260nm 
sample in 25nm < I < 50nm scaling region. The graph 
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FIG. 3. Log-log plot of the second moment of 
height-difference vs I, which shows that for samples 
with thickness 20nm and lOOnm the scaling region is 
4nm < / < 20nm and for the sample with thickness 260nm 
there is a cross over scale /* « 22nm. The upper and lower 
parts of height-height correlation function have the roughness 
exponent 0.29 ± 0.04 and 0.83 ± 0.03, respectively. 

of £ g for the 4nm < I < 2Qnm region is the same as Fig. 5 
but with different slope. In the two scaling regions, the 
£ g 's have a linear dependence on q. This measurement in- 
dicates that the height-fluctuations are not intermittent 
and all of scaling exponents in the stationary state can 
be expressed by £ 2i , i = 1,2 only. 

There is another method to check that the intermit- 
tency is absent in the height fluctuations of V2O5 sur- 
face [12-13]. Let's introduce the generating function 
Z(q,N;l). The generating function is defined through 

N 

Z(g,iV;0 = 5>?, W 

where the normalized measure \ii > is 

\h(xi + I) - h(xi)\ 
^ = — m 1 ( 5 ) 

Y:li\Hxi+i)-h( Xi )\ 

and, N is the number of data \h(xi + I) — h(xi)\. For 
large N, the generating function Z scales as Z(q, N; I) ~ 
N~ Ti - q \ We measure the exponent r{q) = (q — l)D q in 
both scaling region ( 4nm to 20nm and 25nm to 50nm). 
The amplitude of Z depends only on I and q. In Fig. 6 we 
plot the (1 — q)D q vs q for small I. The figure shows that 
the r(q) has a linear dependence on q and the generalized 
fractal dimension D q is independent of q. Its value is 
D q = -Do = 7.34 ± 0.02. Therefore we conclude that the 
height fluctuation is not intermittent. It is necessary to 
note that the exponents £ g and D q are not independent 
in the small region of /. Defining £ g = qH q , it has been 
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FIG. 4. log-log plot of the third moment of 
height-difference vs I, which shows that in the length scale 
I* ~ 22nm there is also a cross over for the scaling exponent 
£3 of the sample 260nm. 
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FIG. 5. The plot of £ q vs q, which has a linear dependence 
on q and shows that the height-difference fluctuations are not 
intermittent. 




FIG. 6. Plot of (1 — q)D q vs q, shows that the fluctuations 
are determined only by dimension Do and the surface are not 
intermittent in the scaling region. 

shown in [20] that H q can be expressed in terms of D q 
as; 



Hq = Hi 



(g -!)(£>,-£>, 



eff) 



q 



(6) 



Our measurement for H„ shows that the H„ 



is inde- 
pendent of q and this means that -D e // = D q = Dq. 
Let us discuss the origin of the -D e // in eq.(6). To 
evaluate the C„ 



we should calculate the summation 
1 \h(xi) — h(xi + l)\ q , where N is the number of 
points over which the average is taken. Normally small 
I is meant by I ~ 1/N. The authors of Ref.[20] assumed 
that, in evaluation of C q , I and TV may be related in a way 



different from l^l/N. That is N - I 



— D, 



( D e jf could 



be considered here as a fractal dimension of the effective 
support of the process). The choice of a particular par- 
tition has no effect on the H q spectrum. However, D e ff 
enters the relation between H q and generalized dimension 
D q as expressed in eq.(6). 



IV. DISCUSSION AND GROWTH MODEL 

In our experiments we obtain a region in which the 
growth exponents of V2O5 are consistent with early- 
scaling regime of forced Kuramoto-Sivashinsky (KS) 
equation in 2+1 dimensions [19]. The forced Kuramoto- 
Sivashinsky has the following form [14-19], 

^- = vV 2 h - k\7 4 h + (A/2)(V/i) 2 + ij{x, t) (7) 

where the v,k and A are the surface tension, surface diffu- 
sion and non-linear factor, respectively. The force 77 is a 
noise term reflecting spatial and temporal fluctuation in 
the incoming flux of material and has a Gaussian distri- 
bution and uncorrelated in space and time. In the limit 
k — and for v > the forced KS equation reduces to 



the Kardar-Parisi-Zhang (KPZ) equation. The KS equa- 
tion with A = is a linear equation and can be exactly 
solved by the standard methods [19,22]. For some Gaus- 
sian noise term, because of the linearity the probability 
distribution function (PDF) of h — h is also Gaussian. 
It is discussed in [22] that for A = and v/k < the 
KS equation generates a mound surface and for v/k > 
gives a self-afine surface, with a roughness exponent ~ 1. 
For A ^ the nonlinear term breaks the symmetry un- 
der transformation h — > —h. Therefore the PDF of h — h 
must be skewed. Our measured value for the PDF of 
h — h shows that the PDF has positive skewness. We will 
report the skewness and kurtosis of the PDF of the V2O5 
films elsewhere [12]. 

Recent simulation of KS equation reveals the presence 
of the early and long scaling regimes [19]. The initial- 
time values for the growth exponent (3 and the roughness 
exponent x are found to be 0.22 — 0.25 and 0.75 — 0.80 re- 
spectively. The long time scaling regime is determined by 
the exponents /? = 0.16 - 0.21 and x = 0-25 - 0.28. The 
scaling exponents are notably less than the exponents 
of KPZ equation [19]. The long-time behavior of the 
KS-cquation has an interesting feature. For long-times 
the height-height correlation function exhibits a bifrac- 
tal structure with two different roughness exponents (see 
Fig.ll in rcf. [19]). For A = 2.0, v = -0.2 and k = 2.0, 
the upper and lower parts of height-height correlation 
function have the roughness exponent 0.27 and 0.71, re- 
spectively. 

Now let us compare our results with the numerical sim- 
ulation of KS equation. Using the numerical results of 
ref. [19] in table II, one can observe that for early-scaling 
regime of the Kuramoto-Sivashinsky equation the expo- 
nents x an( f z satisfy the exponent-identity x + z — 4. 
Using the time dependence of roughness of the samples 
with thickness 20nm , 50nm, lOOnm, 120nm and 150nm 
we find that (3 = 0.29 ± 0.04 (see Fig. 7) and therefore 
X + z ~ 3.7 ± 0.4, which satisfies the exponent-identity 
within the experimental errors. We also observe the sim- 
ilar bifractal structure of the height-height correlation 
function in our samples with thickness 260nm (see Fig.ll 
in [19]). To relate the cross over scale I* to the the sur- 
face morphology, we have measured average distance of 
nearest local maxima d max , and obtain that the d max is 
in the same order of magnitude of Z*. This means that, in 
the stationary state, in average, if the distance between 
the points x\ and X2 lies between the two local maxima 
of h(x), the dynamics is determined by the roughness ex- 
ponent xi- Also for the points that the relative distance 
lies between the next neighbor maximums ( ~ 20 — AOnm 
), the dynamics is given by another exponent X2 which 
is less than the KPZ roughness exponent. Our measure- 
ments show that the roughness exponent xi = 0.29±0.04 
is less than the the KPZ roughness exponent in 2 + 1 di- 
mensions i.e. 0.38. 

In summary we have checked the intermittency of 
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FIG. 7. Log-log plot of the interface width w vs deposi- 
tion time or thickness (in A) t for the samples with thick- 
ness 20nm, 50nm, lOOnm, 120nm and 150nm. The measured 
value for fi is 0.29 ± 0.04. 
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FIG. 8. Log-log plot of the 

second moment of height-difference vs I, for polished Si(100) 
substrate. The measured value for the roughness exponent is 
XS* = 0.47 ±0.03. 




FIG. 9. AFM image of polished Si(100) substrate. 



height fluctuation of V2O5 via two different methods. 
These methods predict that the statistics of height fluc- 
tuations is not intermittent. We measure the scaling ex- 
ponents of height-difference moments, C q (l = \xi — £2!), 
and show that for small I, they behave as ~ \xi — X2^ q - 
The obtained £ 9 's are linear function of q. The observa- 
tion can be explained with the recent numerical simula- 
tion on the basis of nonlinear KS equation. However, one 
should note that for very early stages, when v term domi- 
nates k term, instability arises in x or y direction causing 
the ripple structure with corresponding wave vector in x 
or y direction in KS equation. When the ripple struc- 
ture is formed, there exists slope asymmetry that acti- 
vates the nonlinear effect. Therefore, in the early-stage 
of the process, the ripple structure transforms to sym- 
metric mounds. In our experiments we were not able to 
detect the ripple structure for the samples with < h > 
■C 20nm. It could mean that the polished Si(10Q) sub- 
strate may have an initial roughness which may destroy 
the ripple structure of the film for very early-stages of 
the growth. We measure the roughness exponent of the 
silicon substrate and find that the Xsi = 0.47 ± 0.03 (see 
Figs.8 and 9). 

Finally we note that the measured roughness expo- 
nents x an d P can be higher than the true values because 
of the tip effect [21]. Aue etal. showed that the surface 
fractal dimension (fractal dimension df for 2+1 interface 
is related to % by g?/ = 3 — x ) determined with a scanning 
probe technique will always lead to an underestimate of 
the actual fractal dimension, due to the convolutions of 
tip and surface. Their analysis included tips with differ- 
ent shapes and aspect ratios. For a tip similar to what 
we have used, it is suggested that our true xi an d X2 
should be around ~ 0.75 and ~ 0.25, respectively We 
note that the corrected exponents also satisfy the expo- 
nent identity. 
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